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SUMMARY 

Optimum hat-stiffened compression panel designs are determined using a structural 
synthesis technique. Effects of simplifying assumptions made in the buckling analysis for 
the optimization program are investigated using a more accurate analysis which is a linked 
plate element program. Optimization results for an aluminum panel are compared with 
available results. Optimization results for hat-stiffened graphite -epoxy panels show a 
50-percent weight savings over optimized aluminum panels. Using the structural synthe- 
sis technique, composite panels are shown to possess a variety of proportions at nearly 
constant weight. 


INTRODUCTION 

Optimization of structural members has been a necessary research objective in the 
past (refs. 1 to 12). The need for highly efficient structures in the aerospace industry led 
to the development of fiber-reinforced composite materials. The low density and high 
stiffness of these materials, relative to conventional aerospace metals, indicated that they 
would be highly efficient in compression members, and that the optimum proportions of 
such structures should be investigated. 

Past investigators (refs. 13 to 19) have used various techniques to determine opti- 
mum proportions for structural elements constructed from conventional isotropic metals. 
These techniques, for many reasons, are insufficient for the optimum design of composite 
structural elements. Unlike conventional metals, filamentary reinforced composite mate- 
rials are orthotropic, or even anisotropic. The full determination of the appropriate form 
and optimum proportions of structural elements using these materials requires considera- 
tion of significantly larger numbers of design variables than do isotropic designs. Impos- 
ing undue restrictions on the design process just to effect a simpler solution to the optimi- 
zation problem denies the designer many of the advantages offered by using composite 
materials. Therefore, the purpose of this paper is to present optimized designs for hat- 
stiffened graphite -epoxy compression panels, obtained using an optimization technique 
which can efficiently handle a large number of design variables. 



In this study, a nonlinear mathematical programming technique (ref. 20) called 
AESOP (Automated Engineering and Scientific ^timization Program) was used in conjunc- 
tion with a stiffened-panel mathematical model to determine optimum composite panel 
designs. Optimum stiffened designs were also generated using aluminum properties for 
comparison with the graphite- epoxy results and published aluminum designs. 

The mathematical model of the stiffened panel considered both stability and strength 
of the panel elements. Local stability of the panel elements was described by orthotropic 
plate theory. Overall panel stability (Euler buckling) was described using wide- column 
theory. A simplified maximum strain criteria was used to describe strength limitations 
of the graphite -epoxy material. 

SYMBOLS 

The units used for physical quantities defined in this paper are given in the Later- 
national System (SI) of Units, except where noted. Correlations between this system of 
units and U.S. Customary Units are given in reference 21. 

A total cross sectional area of one pitch of the stiffened panel 

Aj[ area of the ith panel element 

b total panel width 

bj width of the ith panel element 

b^ plate width 

^2 projected depth of panel stiffener 

Dy bending stiffness coefficients 

E Young's modulus 

Ej elastic modulus of the ith panel element 

El effective bending stiffness of the panel 

E^ principal lamina modulus in the direction of loading 
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Ell lamina modulus in the filament direction 

E22 lamina modulus transverse to. the filament direction 

fj percentage of ±45° material in the ith panel element 

G12 lamina inplane shear modulus 

Iq. moment of inertia of the ith panel elements about its centroid 

i,j indices 

L length of the panel 

Nx applied load per unit width of the panel 

NxEuier Euler buckling load per unit width of the panel 

P„ applied load on the ith panel element 

local buckling load of the ith panel element 
P total load on the panel per unit pitch 

t plate thickness 

tj thickness of the ith panel element 

t^ lamina thickness 

W panel mass 

y distance of the neutral axis from the reference axis 

Yj distance of centroid of the ith panel element from the reference axis 

e strain in the panel 

€ 3 , allowable strain 
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V 


Poisson’s ratio 


1/12 lamina transverse Poisson’s ratio 

p material density 

p mass per unit area 

(f) performance function 

CT„ stress applied to the ith panel element 

local buckling stress 

Subscript: 

i ith panel element 

Superscripts; 

L lower bound 

H upper bound 


ANALYSES 

The general optimization cycle used is depicted in figure 1. The synthesis model 
consists of two parts, namely, the mathematical model and the optimizer. The mathe- 
matical model describes the strength and stabili;y idealization of the structure to be opti- 
mized and a structural analysis of the mathemat ^.l model determines the design varia- 
bles. The performance function is evaluated for hese design variables and checked for 
optimality. The design variables are incremented, in accordance with the optimization 
Scheme, subject to various constraints and the process is repeated until an optimum 
design has been generated. Details of the optimization scheme are given in appendix A. 

The problem considered herein is the optimization of a hat- stiffened composite panel 
under uniaxial compression. The loaded edges of the panel are simply-supported and the 
unloaded edges are free. Figure 2 shows a schematic drawing of such a panel. Proper- 
ties of the graphite-epoxy and aluminum material used in this study are given in table 1. 
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Mathematical Model 


Basic assumptions .- The following basic assumptions are noted: 

(1) Panel elements for the local buckling analysis are orthotropic constant-thickness 
plates, simply supported on all four edges. 

(2) The panel is assumed to behave like a wide column for the overall or Euler buck- 
ling analysis. 

(3) Twisting (or torsional) failure modes of the stiffeners are not considered. 

(4) Two kinds of orthotropic lamina are used, namely, 0° lamina and a lamina with 
the extensional properties of a ±45° laminate. 

(5) Bending-twisting coupling effects in the composite material are not considered. 

(6) Each panel element is assumed to have no more than three layers, stacked in a 
balanced and midplane symmetric manner. Either 0° or ±45° material is permitted as the 
outside laminae. 

(7) Allowable strain in compression for any panel element is conservatively set equal 
to the yield strain of the 0° composite material irrespective of the percentage of ±45° 
material contained in the panel element. 

Design variables .- Under the assumption of wide- column behavior, only one pitch of 
the stiffener spacing is required for analysis. Figure 3 shows a representative cross sec- 
tion of the idealized panel and details the design variables used. For the hat-stiffened 
configuration the design variables are: 

bj width of panel element 

tj thickness of panel element 

fi percentage of ±45° material in panel element where i = 1, 2 , 3, 4 

Performance function .- The performance function used in this analysis is a weight 
parameter defined as mass per unit width of the panel 

. = PA 

(bi + 2b4) (1) 

Constraints . - hi order to complete the definition of the problem, geometric and 
strength constraints must be specified. Those selected for this study are: 
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(1) Local buckling load of each panel element shall be greater than or equal to the 
applied load, 




(2) Euler buckling load of the total panel shall be greater than or equal to the applied 
loading, 

N„ s Nx (3) 

^Euler ^ 

(3) Applied strain of the total panel shall be less than or equal to the allowable strain. 


(4) Stiffener spacing shall be greater than or equal to bg (see fig. 3), 


bj + 2b^ >bg 


(5) The value of the design variables shall be limited to a region of practical interest. 


bf* sb. sb? 

1 1 1 


t^ < t ^ t^ 
1 *^1 


fL < f < fH 
i i i 


For this study these limits are: 


0.762 s bj ^ 25.4 cm 


0.00254 <tj £ 1.27 cm 


0 f j ^ 100 percent 


The values of the load and strain parameters in equations (2) to (4) are determined through 
a simplified buckling analysis discussed in detail in appendix B. 
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RESULTS AND DISCUSSION 


Using the synthesis model described above, optimum designs for graphite -epoxy 
and aluminum panels are generated at various load levels. Results are summarized in 
figure 4, which is a standard weight-strength plot of the strength parameter Nx/L and 
the weight parameter W/bL^. Results for the composite and aluminum panels are shown 
as solid lines. The broken lines are the material strength limits for the aluminum and 
the graphite -epoxy material. 


Aluminum Panels 

Results from reference 2 for an aluminum panel are presented in figure 4 for com- 
parison. As shown, the present results for an aluminum panel show a slight advantage 
over the results of reference 2. In reference 2 the flange width, cap width, and stiffener 
height are arbitrarily set to be a percentage of the width bj of the stiffener to effect an 
easier solution. No such simplifying assumptions with respect to the cross section are 
used in the present study, hence the optimum panel designs show a slight advantage over 
those of reference 2. 

Also, reference 2 uses the optimization condition that local buckling in each panel ele- 
ment be coincident with Euler buckling of the whole panel. In the present analysis no such 
condition is imposed. However, results of the optimization process show that, in fact, 
for the optimum aluminum design, local buckling of each element and Euler buckling of 
the entire panel should occur simultaneously. 

Graphite-Epoxy Panels 

An attempt was made to compare results obtained for graphite -epoxy panels with 
similarly optimized composite panels in the literature, however, no published results 
could be found. Therefore, to assess the accuracy of the procedures used, a study is 
made of the effects of the simplifying assumptions employed in the present buckling analy- 
sis. This is done by determining the buckling load for a set of optimized panels using 
BUCLASP-2 (refs. 22 and 23), a computer program which is devoid of all such assump- 
tions, with the exception that bending-twisting coupling is ignored. The BUCLASP-2 buck- 
ling loads are then plotted against the optimum value of the weight parameter. These 
points, shown as circles in figure 4, show good agreement with the present results. It is 
also noted that for the composite panel, results obtained using BUCLASP-2 show that local 
buckling and Euler buckling of the optimized panel occur at very nearly the same load. 

The results in figure 4 show that optimized graphite-epoxy panel designs weigh 
approximately one-half as much as optimized aluminum designs over a wide loading range. 
Thus, such panels have a significant potential for lightweight aircraft or spacecraft com- 
pression structures. 
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An important result of the present study is illustrated in figure 5 which shows two 
composite panel cross sections that weigh approximately the same at both a high and a low 
load. The design variables for a range of loadings that show this phenomenon are detailed 
in tables 2 to 5 and differ markedly. This result has fortuitous ramifications when manu- 
facturing constraints are considered, indicating that some design variables could be fixed 
for manufacturing reasons without adversely affecting the efficiency of the panel. 

However, one should make sure in such cases that the extreme designs do not come 
about because of assumptions made during the buckling analysis and, in reality, fail in a 
neglected mode. For example, figure 6 shows the buckling mode shapes for the two highly 
loaded panels in figure 5. These mode shapes were obtained from BUCLASP-2. The 
panel in figure 6(b) fails in a local buckling mode and all the panel members behave as 
simply supported elements, which is in accord with the simplifying assumptions. The 
stiffener in figure 6(a) is very deep, compared to its width, and fails in a twisting, or 
torsional mode, 10 percent below the optimum design value. This mode of failure was 
neglected in the mathematical-model stability analysis. Alternate designs which are very 
similar to that of figure 6(a) but do not exhibit the torsional mode of failure can be found 
in tables 2 to 5. 

The panel shown in figure 5(b) has 100 percent ±45° material in the skin and in the 
stiffener webs. Also, the optimized thickness of the webs and skin is very small compared 
to the thickness of the stiffener cap and flanges (see tables 2 to 7). This suggests that 
most of the load is carried by the 0° material, which results in a more efficient panel. It 
is significant that, in general, the optimization procedure found designs with 100 percent 
0° material in the flanges and stiffener caps and 100 percent ±45° material in the webs. 
This suggests a method of reinforcing conventional metal hat-stiffened panels with com- 
posites by bonding unidirectional material to the flange and stiffener cap. Li practice this 
is indeed the case if the composite material is compatible with the metal in the panel. It 
should be noted in figure 5 that the amount of ±45° material in the skin ranges from 0 to 
100 percent. Only the two extreme cases, however, are shown in the figure. 

Practical Constraints 

Results presented in figure 4 and tables 2 to 8 give the optimized values of the dif- 
ferent design variables for the composite and the aluminum panels; examination of these 
results reveals that they are not very practical from a manufacturing point of view. In 
the case of the aluminum panels it is more practical to have a constant stiffener thickness 
(i.e., tg = tg = t^ (see fig. 3)). The effect of imposing such a constraint in the optimiza- 
tion procedure on the aluminum panel is shown in figure 7. It can be seen that an approxi- 
mate 10-percent weight penalty is paid in using such a constraint. 

In the case of the composite panel it is not practicable to make a hat stiffener with 
only 0° material in the stiffener cap. A more practical design, shown in figure 8, 
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eliminates the problem of bonding the 0° cap material to the ±45° material in the webs. 
The 0° material for the cap is encapsulated by the ±45° material which makes up the stif- 
fener webs. The effect of placing such a constraint on the optimum panel is shown in fig- 
ure 9 and, as can be seen, only a very nominal penalty is incurred. If practical con- 
straints are included in the optimization process, the composite material appears to give 
more flexibility to the design process with little penalty. Values of the optimized design 
variables for the composite panel with practical constraints are presented in table 9. 
Specifying the orientation of the material comprising the various plate elements, as shown 
in figure 8, eliminates fj as a design variable and reduces the total number of variables 
from 12 to 8. 


Variable Sensitivity 

An apparent insensitivity of the weight parameter W/bL^ to design-variable change 
was investigated in detail to define limits for certain design variables. Results of studies 
of the variables bj and b 3 are shown in figures 10 and 11, respectively. These fig- 
ures are obtained for the fixed value of N^/L equal to 2069 kPa. 

The plot of W/bL^ and bj in figure 10 is obtained by specifying bj during an 
optimization cycle while the other design variables remain unspecified. Since previous 
optimization results always give the same values for f 2 , f 3 , and f^, these values were 
also held fixed, leaving fj free to vary. A similar approach was used for the results in 
figure 11. Figures 10 and 11 show that the weight parameter is relatively insensitive to 
independent changes in either bj or b 3 over a wide range (shown by the hatched area). 
This demonstrates that considerable freedom exists for a designer to specify practical 
dimensional constraints without incurring any weight-parameter penalty. Optimum con- 
struction for the panel skins was also found to be a function of bj and b 3 . The variable 
f j varied from 0 percent for the narrow stiffener design to 100 percent for the wide 
stiffener. 


CONCLUDING REMARKS 

A mathematical model including both strength and stability effects, for hat-stiffened 
compression panels has been developed. This model was combined with a nonlinear math- 
ematical programming technique and used to generate optimum designs for both graphite- 
epoxy and aluminum panels. Selected optimum designs were analyzed with a more com- 
plete stability analysis which was devoid of assumptions used in the synthesis model. 

Based on the analytical results presented herein, the following conclusions can be made: 

1. Optimized graphite-epoxy panels weigh approximately one-half as much as opti- 
mized aluminum panels over a wide loading range. 
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2. Composite panel designs can be evolved via the optimization procedure which 
weigh essentially the same, but have geometric proportions which differ markedly. 

3. Introduction of practical constraints as design variables can in some cases result 
in panels which are easier to manufacture and which will only weigh 10 percent more than 
optimum. 

4. In order to ensure that complex failure modes will not cause premature failure, 
optimized designs should be reanalyzed with methods which are devoid of the simplifying 
assumptions incorporated in the synthesis model. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 16, 1974. 
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APPENDIX A 


DESIGN SYNTHESIS PROGRAM 

This appendix is devoted to discussing the design synthesis program developed for 
the purpose of optimizing the weight parameter of a composite hat-stiffened panel under 
uniaxial compression. 

Figure 12 shows a schematic flow chart of the optimization computer program AESOP 
(Automated Engineering and Scientific ^timization Program). (See ref. 20.) There are 
nine search algorithms available in AESOP, for the purpose of seeking a minimum solu- 
tion. One or a combination of search algorithms can be used for this purpose. For the 
present work a combination of adaptive creep search and pattern search was found to be 
most efficient and effective in terms of computation and, hence, was used for the optimi- 
zation work. Resulting typical CPU times for a single panel-optimization run on the 
CDC 6600 computer were about 10 seconds. 

The AESOP computer program requires a user-defined subroutine describing the 
problem to be minimized. The hatched area in figure 12 represents subroutine PANEL 
developed for the present problem. A listing of subroutine PANEL is given at the end of 
this appendix, fi should be noted that subroutine PANEL has been developed only for the 
composite hat-stiffened panel. In order to use this program for any other kind of panel it 
would have to be modified. 

Subroutine PANEL employs the user-defined subroutines SIGEULR and SIGLOCL 
which compute the Euler buckling stress and the local buckling stress, respectively. 
Listings of these subroutines are also given with the listing of PANEL. 
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1800129 


EX(NN)=SUM5 

1800130 


0CNST( NN) = S0PT( SIIM1*SUM2 ) «-SUM3*-2 . *S()M4 

1800131 


TF(NN.EQ.2) fi(2)=R(6) 

1800132 


TF (NN.NF.3) GO Tn lo 

1800133 


TT(3)=TT (3) *-TT(2) 

1800134 

10 

IFINN.NE.A) GO TO 1 

1800135 


TT(4)=TT(4|fTT(l).TT(?( 

1800136 


« (4) =2.*B( 4 ) 

1 800137 

I 

P TGMA{ NN ) = AK ( NM ) «<p '**?*( ( NN) /R ( NN) ) «*?<=Cr''"$T (NN ) /o . 

1800138 


R ( 2) =RAA 

13001 39 


TF (NN.MF.3) TO 20 

1800140 


TT(i)=TT(3)-TT(2) 

18U0141 

?0 

TF(N'‘.Nt.4) GO TO 2 

1800142 


tt(41=tt(4)-TT( 1 )-tt(2 ) 

1300143 


R(4) = B(4 M7. 

1800144 

7 

pFTUbn 

1800145 



FA'D 

1800146 
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APPENDIX A - Concluded 


c T \if-- c; inrui ^ 


1300147 



1300148 

'■rMV 1N/HTP/R ( 10 ) , TT ( 5 ) , FX( 6) , fXL .f' ' , ';iGM4< lU) . Ak tA 


1300149 

M i) = A( 1 |*TT( I ) 


1300150 

* AA>tTT ( ? 1 *HX ( 2 1 / r X ( 1 1 


1800151 

A n 1 = Ft ( 1 ) * ( T T ( 3 ) + T T ( ? ) ) * F X ( n / t X ( u 


1800152 

A(/+1=?.*p(4I*(TT(1M-TT(f, )fTT(2n*EX(Al/FX{l) 


1800153 

A (5 1 =A( n f A ( 2 ) f A ( 3 n-A ( 4) 


1800154 

VAAP = (A(2l*H(2)/2.+A(3)*R(2)*.A(4l#(’'’’(l>+TT(4)»TTi2))/2. 

) /A ( 4 1 

18001 55 

A(7)=A(3)/(TT(3)fTT(2)l 


1800156 

A (i! 1 =4( 4 1/ ( 2 .*‘FT ( 4 1 +TT (? IfTT ( 1 ) ) 


18001 57 

’T ( 4) =tt (^)*FX(2)AFX(1I 


1300158 

Xkri-A(ll>'T'F(ll*>"3/12.*-A(cJI*(TTtl)HT(2)«-TT(4ll*#3/t,.H5(7)#(TT(3) 

1800159 

1 tTT(3||**3/12.+‘FT(5)*(‘>(2)-(''’T(l)-TT(2)-TT(^l)/2.)#w 


1800160 

2 3/6. fA( 1 ) pyP A9**2 F A ( ? ) *( Y RAR-3( 2 ) / 2. ) #*2 *-A (3 )*4YF.Ak-8 1 ^ ) 

)**2<-A( 4)* 

1800161 

3(YRA9-(TTm+TT(4)tTT(2))/2.)**2 


1800162 

^ Tr.vA ( 6 ) =F X ( ]. ) - P ^ wnr / ( A AP F C ) 


18001 63 

' J P M 


1800164 

r N’n 


1800165 
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APPENDIX B 


BUCKLING ANALYSIS 

The buckling analysis used is based on the variables for the cross section shown in 
figure 3. A similar analysis for the case presented in figure 8 can be performed by rede- 
fining the corresponding variables. 


Load in Each Panel Element 


Assuming Nx is the load intensity per unit width, a 

i 

panel element, and P is the total load per stiffener spacing. 


is the axial stress in each 
then. 


P = Nx • (bi + 2b4) 


(Bl) 


and 


P = 


4 


Z 

i=l 



For compatibility, the strain in each panel element has to be equal. Hence 


e 



(i = 1, 2, 3, 4) 


(B2) 


(B3) 


where 


and E^ is the principal lamina modulus calculated from the properties in table 1 through 
a standard lamina transformation formula (ref. 24). The Poisson effect neglected by this 
approach was found to be less than 0.5 percent. From equations (B2) and (B3) 

4 

P = Z 

i=l 

Solving this equation for the strain e and using equation (Bl) 
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APPENDIX B — Continued 


e = 


Nx(bi + 2b4) 

“4 

Tj ®i-^i 

i=l 


Finally the load Pa^ in each panel member is then given by 


NybEiAi 

I’a, = <’a A, = 


Z ElAi 

i=l 


(B4) 


(B5) 


Local Buckling 

Each panel member is assumed to be orthotropic and simply supported on all four 
edges. Hence, from orthotropic plate theory (ref. 25) 


2 

(\/^ii^22 ^12 + ^^ee) 

b^t ^ ' 

and 

P^ = cr^b^t 


where, 

cr^ local buckling stress 

b^ width of the plate 

t thickness of the plate 

Djj bending stiffness coefficients 

P^ local buckling load 

Here for the sake of simplicity the subscript i has been omitted. However, this equation 
applies to each panel element. 
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APPENDIX B — Continued 


Euler Buckling 

Only one pitch of stiffener spacing need be considered for the purpose of the Euler 
buckling analysis. Since each panel element can have a different percentage of ±45° lami- 
nates, each panel element will have a different value of Young's modulus. Ih order to find 
the Euler buckling load the equivalent area approach is used to find the effective Young's 
modulus and effective area of each panel element. Then 


N- 


^Euler 


7r^(EI) 


bL 


(B7) 


where El is the effective stiffness. 

Assuming y^ is the distance of the center of gravity of the ith panel element from 
the reference axis (shown in fig. 3), E. is the principal lamina modulus for the ith panel 
element, A^^ is the area of the ith panel element, and y is the distance of the neutral 
axis from the reference axis, then, 

4 

Z EiAjyi 

y = ^ (B8) 

Z EiAi 

i=l 

From figure 3 it can be seen that 



A4 = 2b4(ti + t4) 


Hence, 
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APPENDIX B - Concluded 




y = 




E4A^(tj + t^) 


4 

z 

i=l 


EiAi 


(BID 


The effective El of the cross section about the neutral axis can be written as 
4 4 

El= I Ejlo. + Z EiAiYi^ (B12) 

i=l ^ i=l 

where 

lo moment of inertia of the area of the ith panel element about its centroid 

1 

Yj distance of the centroid of the ith panel element from the neutral axis of the 

cross section 

Hence, 

El = ± Ei(b, t 2b4)ti3 + i Ejtj 

H ^ 2b4)tjy^ 

+ E2A2(y-ib2-^) + EjAjS - bjj^ 

+ 2E4b4t4|? - I fti + tj (BIS) 


2br 


2bg - - t^ 


^^2 ■ ^4) 
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APPENDIX C 


BUCLASP-2 ASSUMPTIONS AND MODEL 

Effects of the simplifying assumptions in the buckling analysis were studied using 
the more accurate linked plate analysis program BUCLASP-2. This appendix is devoted 
to a discussion of some of the capabilities of BUCLASP-2 (a computer program for the 
instability analysis of biaxially loaded composite panels) as it pertains to the buckling anal- 
ysis of the composite panels considered in the present work. This computer program 
(refs. 22 and 23) is operational on the CDC 6600 computer. Some of the basic assump- 
tions made in the analysis of BUCLASP-2 are as follows: 

(1) Panel elements are orthotropic and have balanced laminates. 

(2) The material is linearly elastic. 

(3) Thin-plate theory is employed. 

(4) Effects of prebuckling deformations are ignored. 

(5) Edges normal to the longitudinal direction are assumed to be simply supported. 

Support conditions at other boundaries are arbitrary. With the above assumptions 
an "exact" analysis of the whole panel is made. This analysis results in the prediction of 
Euler buckling modes, local buckling modes, or coupled Euler and local modes. 

The user of BUCLASP-2 has to define the mathematical model of the panel under 
consideration. This mathematical model consists of three substructures, namely, the 
start substructure, end substructure, and the repeat substructure. Figure 13 shows the 
cross sections of the three substructures for the panel studied in this investigation. The 
optimized panels studied using BUCLASP-2 in this work were seven stiffener spacings 
wide. The results after using AESOP define the cross-sectional dimensions of the panel. 
These dimensions are then used to find the buckling load using BUCLASP-2. 
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TABLE 1.- MATERIAL PROPERTIES 


(a) SI Units 


Graphite epoxy 

Aluminum 

Ell = 138 GPa 
E22 = 8.962 GPa 
G12 = 4.48 GPa 
1^12 - 0.304 
p = 1522 kg/m^ 

= 0.00575 

E = 68.9 GPa 
u = 0.300 
p = 2768 kg/m^ 
= 0.00575 


(b) U.S. Customary Units 


Graphite epoxy 

Aluminum 

Ell = 2 X lo”^ psi 
E22 = 1.3 X 10® psi 
G12 = 6.5 X 10® psi 
1^12 = 0.304 
p = 0.055 Ib/in® 

= 0.00575 

E = lo'^ psi 
u = 0.300 
p = 0.100 Ib/in® 
ea = 0.00575 






































TABLE 4.- OPTIMIZED DESIGN VARIABLES FOR GRAPHITE -EPOXY PANEL 


fl = fg = 14 = 0; fg = 100; 


L = 127.0 cm 



(a) SI Units 


B 

JV_ 

bL^’ 

kg/m3 

bi, 

cm 

B 

H 

b4, 

cm 

cm 

tg, 

cm 

*3» 

cm 

t4> 

cm 

— 


143.8 X 10"^ 


6.401 

3.099 

1.651 

0.549 


1.082 


5516 

124.4 


6.502 

1.930 

1.854 

.500 

.102 

1.270 


3448 

91.2 

3.378 

6.350 

1.270 

2.565 

.447 


1.270 

.025 

2758 

83.5 

3.378 

5.563 

3.175 

2.743 

.287 


.559 

.178 


69.1 

2.972 

4.826 

2.540 

1.651 

.203 


.356 

.025 

1379 

55.9 

2.718 

4.724 

.762 

1.600 

.206 


1.245 

.036 



2.362 

4.013 

.762 

1.448 

.147 

.043 

.737 

.036 

345 

28.5 


3.835 

.762 

1.168 

.107 

.041 

.267 

.025 


(b) U.S. Customary Units 


Q 


99 




ti, 

in. 


in. 

t4, 

in. 

1000 


1.340 

2.520 


0.650 

0.216 

0.040 

0.426 

0.010 

800 



2.560 


.730 

.197 

.040 

.500 

.027 

500 


1.330 


.500 


.176 

.038 

.500 

.010 

400 


1.330 

2.190 



.113 

.033 

.220 

.070 

300 


1.170 



.650 

.080 

.025 

.140 

.010 

200 



1.860 

.300 

.630 

.081 

.023 

.490 

.014 

100 


.930 

1.580 

.300 

.570 

.058 

.017 

.290 

.014 



.750 

1.510 



.300 

.460 

.042 

.016 

.105 

.010 
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TABLE 5.- OPTIMIZED DESIGN VARIABLES FOR GRAPHITE-EPOXY PANEL 


[fj = fg = 100; fg = = 0; L = 76.2 cm (30 in.j 

(a) SI Units 


Nx 
L ’ 
kPa 

W 

bL2’ 

kg/m3 

bl, 

cm 

^2^ 

cm 

b3. 

cm 

b4, 

cm 

tl, 

cm 

t2> 

cm 

cm 

4, 

cm 

3448 

92.6 X 10"^ 

2.718 

4.267 

1.397 

2.337 



1.270 

BSi 

2758 

81.6 

4.318 

4.115 

3.861 



.076 

.434 


2069 

66.4 

4.140 

4.140 

1.626 

2.159 

.051 

.081 

.549 

■9 

1379 

57.5 

2.350 

3.277 

2.108 

2.388 

.038 

.051 

.345 

.152 

690 

39.0 

2.197 

2.946 

.762 

1.712 

.036 

.041 

.470 

.091 

345 

27.4 

1.999 

2.565 

1.651 

1.321 

.025 

.033 

.127 

.058 


(b) U.S. Customary Units 


m 

W 

bL2’ 

lb/in3 

BH 

■1 

H 

BH 

■1 

m 

4’ 

in. 

tg, 

in. 

in. 

in. 

500 

3.350 X 10'^ 

1.070 

1.680 

0.550 

0.920 

0.021 

0.028 

0.500 

0.068 

400 

2.950 

1.700 

1.620 

1.520 

1.000 

.033 

.030 

.171 

.064 

300 

2.400 

1.630 

1.630 

.640 

.850 

.02 

.032 

.216 

.060 

200 

2.080 

.925 

1.290 

.830 

.940 

.015 

.020 

.136 

.060 

100 

1.410 

.865 

1.160 

.300 

.674 

.014 

.016 

.185 

.036 

50 

.990 

.787 

1.010 

.650 

.520 

.010 

.013 

.050 

.023 
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TABLE 8.- OPTIMIZED DESIGN VARIABLES FOR ALL-ALUMINUM PANEL 

WITH L = 76.2 cm (30 in.) 


(a) SI Units 


Nx 
L ’ 
kPa 

W 

bL^’ 

kg/m3 

bl, 

cm 








3448 

230.3 X 10"^ 

8.687 


7.518 

4.775 




0.025 

2758 

185.8 

7.239 

4.242 

6.782 

3.924 

.284 


.287 

.025 

2069 

153.7 

6.142 

4.191 

5.525 


.231 

.150 


.025 

1379 

123.0 

4.572 

3.744 

4.432 

2.642 

.157 

.122 


.025 

690 

87.1 

4.343 

3.226 

3.912 

2.591 

.124 

.089 

.117 

.025 

345 

61.1 

3.061 

2.787 

2.946 





.025 


(b) U.S. Customary Units 


Nx 

L ’ 
psi 

CO 

bD 

in. 


bS’ 

in. 

■ 





500 

8.330 X 10"^ 

3.420 


2.960 

■ 

0.139 

0.058 

0.199 

0.010 

400 

6.720 

2.850 


2.670 

IQ 

.112 

.063 

.113 

.010 

300 

5.560 

2.418 

mm 

2.175 

IQ 

.091 

.059 

.082 

.010 

200 

4.450 

1.800 

mm 

1.745 

llil 

.062 

.048 

.069 

.010 

100 

3.150 

1.710 

IQ 

1.540 

IQ 

.049 

.035 

.046 

.010 

50 

2.210 

1 

1.205 

IQ 

1.160 

IQ 

.029 

.026 

.029 

.010 
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TABLE 9.- OPTIMIZED DESIGN VARIABLES FOR GRAPHITE- EPOXY PANEL 
WITH PRACTICAL CONSTRAINTS 

[^1 "" ^3 " ^4 = ^2 " ^ 

(a) SI Units 


Nx 
L ’ 
kPa 

W 

bL2’ 
kg/m 3 

— 

bi, 

cm 

^2’ 

cm 

bg, 

cm 

b4, 

cm 

^1’ 

cm 

h’ 

cm 

t3, 

cm 

^4’ 

cm 

3448 

92.4 X 10"^ 

3.355 

4.044 

2.289 

3.970 

0.065 

0.072 

0.569 

0.276 

2758 

82.9 

3.091 

3.833 

1.869 

3.241 

.060 

.068 

.632 

.182 

2069 

72.1 

3.719 

3.584 

1.687 

3.155 

.067 

.060 

.620 

.170 

1379 

59.1 

1.918 

3.119 

.869 

2.360 

.030 

.044 

.635 

.153 

690 

40.4 

1.806 

2.786 

.772 

2.090 

.025 

.036 

.414 

.088 

345 

28.7 

2,065 

2.441 

.780 

1.966 

.025 

.037 

.270 

.057 
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Figure 1.- Optimization cycle. 

















100% ±45 
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Figure 5.- Examples of optimized graphite- epoxy sections. 








Design Variables 

t. 4 
1 


Total Number - 8 



±45° material 



material 



Figure 8.- Representative cross section of the panel with practical constraints. 
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strength parameter, N^L, kPa 



Strength parameter, N^L, psi 


Figure 9.- Effect of practical constraints on graphite-epoay panel 
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Figure 13.- Mathematical model for BUCLASP-2. 
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